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1. Introduction 

The experimental accuracy reached by the elec- 
troweak precision observables in particle acceler- 
ators allows to test the Standard Model (SM) 
at the quantum level within an accuracy better 
than one per mil. Although precision observ- 
ables at LEP have played an important role in 
this respect, a new generation of experiments at 
the Tevatron (II), LHC and a future Linear Col- 
lider (LC) will be necessary to find out the ulti- 
mate nature of the spontaneous symmetry break- 
ing mechanism. When extending the Higgs sec- 
tor, we mainly focus on two-Higgs-doublet mod- 
els (2HDM), and most particularly on the Type 
II ones -like that of the the Minimal Supersym- 
metric Standard Model (MSSM) [ |. The fact 
that the Higgs bosons of the MSSM (and their 
SUSY counterparts) couple most preferentially to 
heavy quarks suggests that the physics of the top 
and bottom quarks is the natural arena where to 
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try to get a hint of the existence of SUSY Higgs 
particles. In this paper we report on the produc- 
tion of a heavy charged Higgs boson in association 
with the top and bottom quarks at the Tevatron 
and LHC energies within the MSSM. The mass of 
the charged Higgs is assumed M H ± > 200 GeV, 
and so it cannot come from the decay of a top 
quark. While the detection of a charged Higgs 
boson would still leave a lot of questions unan- 
swered, it would immediately offer (in contrast to 
the detection of a neutral one) indisputable evi- 
dence of physics beyond the SM. The basic pro- 
duction processes under consideration are 

PPipp) H+tb + X Tevatron (LHC) . (1) 

There are many other processes for charged Higgs 
production, which have been studied in the liter- 
ature [|^|, such as: i) one- loop gluon fusion (g g — > 
H + H~); ii) tree-level pair production in bottom 
quark scattering (bb — > H + H~), iii) associated 
production with H /± : bb — > W H T (tree-level) , 
g g — » W ± H T (one- loop) , iv) Drell-Yan type pro- 
duction (qq — > V — > H + H~). However, in all 
these cases the production cross-section is smaller 
than for (g) if we restrict ourselves to high values 
of tan/3 (> 20). Although process (Q) has been 
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studied previously [§, 0, the inclusion of the lead- 
ing quantum SUSY effects was missing and turns 
out to be very important, as first hinted in [ ^| 
and further elaborated in [ ||. 

2. Leading order Cross-section in QCD 

At the parton level, the reaction ([!]) proceeds 
through three channels: i) (^-annihilation for 
light quarks 

qq -> H+ib, (2) 

where q — u,d (the s contribution can be safely 
neglected), a channel only relevant for the Teva- 
tron; ii) gg-fusion 

gg^H+ib, (3) 

which is dominant at the LHC, but it can also 
be important at the Tevatron for increasing H+ 
masses ; and finally there is the iii) bottom-gluon 
2-body channel 

bg -» H+i . (4) 

We will compute the cross-section for the charged 
Higgs boson production process (Q) at the leading 
order (LO) in QCD, namely at 0{a 2 s ). However, 
the QCD corrections at the next-to-leading order 
(NLO) or could be important. Although 

a dedicated calculation of these higher order ef- 
fects for the process under consideration is not 
available, let us mention that the corresponding 
calculation for the related process in the Standard 
Model, pp(pp) — > Hti+X, is available and shows 
that the NLO effects lead to a non-trivial QCD 
A-factor which is typically smaller than one at 
the Tevatron, and up to 1.4 at the LHC [ [|. On 
the other hand, partial calculations of the QCD 
corrections to the process (|l|) at the NLO have 
appeared. For example, in Ref. [ Q the stan- 
dard NLO QCD corrections to the subprocess (|]) 
at the LHC are computed, obtaining a large K- 
factor between ~ 1.6 and ~ 1.8 for tan/3 > 20. 
Taking into account these considerations, a large 
A-factor for the full process ([!]) coming from the 
standard QCD corrections (e.g. K^ cr> ~ 1.5) is 
not ruled out. However, in the absence of a com- 
plete calculation of these QCD effects, in what 
follows we present our results at the leading order 



and will parametrize our ignorance of the com- 
plete NLO effects in terms of a A"-factor, which 
will be assumed different from one only for the 
Tevatron, where the size of the cross-section is 
hardly the necessary for a viable study. On the 
other hand, for the LHC the cross-section, being 
typically three orders of magnitude larger, can af- 
ford a QCD A-factor of one. 

Once a PDF for 6-quarks is used, there is some 
amount of overlap between bg- and gg-initiated 
amplitudes, which has to be removed [ The 
overlap arises because the ^-density in the bg am- 
plitude receives contributions from gluon splitting 
which was already counted in the gg amplitude, 
so we have to avoid double counting by the sub- 
tracting of the gluon splitting term. The net par- 
tonic cross-section from the bg- and gg-initiated 
subprocesses is 

<r(h + 99 -» H+i + X) net = 
a{gb H+i) + a(gg -> H+ib) 

-a(g -^bb®gb^ H+i). (5) 

Furthermore, for the study of the various differ- 
ential distributions one has to properly combine 
the 2—>2 and 2 — > 3 processes mentioned above 
in order to reproduce not only the total cross- 
section but also the correct event kinematics. The 
point is that we know the total amount of double 
counting but not a priori which part of this value 
should be subtracted from the H+i process and 
which part from the H+ib one. We apply here the 
method proposed in [ |j| for the analogous process 
of the single top quark production. According 
to this method we use the cut on the transverse 
momenta of the 6-quark associated with charged 
Higgs boson production to separate and recom- 
bine gg- and 65-initiated processes. 

3. Signal versus background 

Concerning the signature, we consider the dom- 
inant tibb final state for the combined signal pro- 
cess (Q). We focus here on the triple 6-tagging 
case and consider the situation where one top de- 
cays hadronically and the other leptonically (in- 
cluding only electron and muon decay channels) 
in order to reduce the combinatorics when both 
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a(qq -> ttbb) 6.62 fb 0.266 pb 

a{gg -> 0.676 fb 6.00 pb 

cr(flr6 -> ttb) 1.22 fb 4.33 pb 

Subtr. term 0.72 fb 2.1 pb 

Table 1 

The main background processes to the signal ([!]) 
at the Tevatron (2nd column) and the LHC (3rd 
column) under cuts explained in the text. The 
corresponding subtraction term is also shown. 



P M 2 rrig A t A b 

A -1000 200 1000 1000 1000 500 500 

B -200 200 1000 500 500 500 500 

C 200 200 1000 500 500 -500 500 

D 1000 200 1000 1000 1000 -500 500 

Table 2 

Sets of MSSM parameters used in the computa- 
tion of the SUSY corrections to process (|l|). All 
masses and trilinear couplings in GeV. 



top quarks are reconstructed. The branching ra- 
tio of tibb -> bbbbt^vqq 1 is 2/9 x 2/3 x 2 = 8/27. 

In order to decide whether a charged Higgs 
boson cross-section leads to a detectable signal, 
we have to compute the background rate. The 
main QCD backgrounds leading to the same tibb 
signature and their respective cross-sections are 
shown in Table [l]. For reconstructing the tib(b) 
final state, the following steps are followed: i) We 
reconstruct the W-boson mass from lepton and 
neutrino momenta: = (pg+p^) 2 . The basic 

cuts for the lepton (electron or muon) has been 
chosen as follows: 

p e T > 15 GeV, M < 2.5, p% tss > 15 GeV ; (6) 

ii) We reconstruct the mass of the second W- 
boson (M^j). The following basic cuts for the 
jets were chosen for the Tevatron (LHC): 

p>i h > 20 (30) GeV, \ Vj \ < 3, \ Vb \ < 2 . (7) 

Then we form and keep all m tl — Mwu an d 
m t2 = Mw2b combinations for the first and sec- 
ond top-quarks; iii) In the final step we form the 
X function 

X 2 = [Mffi - M w ) 2 + (M^ 2 - M w ) 2 

+ (m tl - m t ) 2 + (m t2 - m, t ) 2 (8) 

for all combinations of 6-jets, jets, lepton and 
neutrino and choose the combination giving the 
smallest (best) value of the \ function. After the 
reconstruction of the ttbb state one should recon- 
struct the charged Higgs boson mass for the signal 
and the continuous tb mass for the background. 
We assume that the 6-jet with the highest pt in 
tbt(b) signature comes from the H + decay. Af- 
ter all cuts are set up we obtained the signal 



and background efficiencies, S/B ratio and signal 
significance S/ \/B— see Tables 4,5 in the second 
paper of [ || for details. They will be used in the 
next section after including the SUSY corrections. 

For the calculation we use m t = 175 GeV, 
m b = 4.6 GeV and the CTEQ4L set of PDFs [ 
[To| . Here m t ,m& refer to the quark pole masses. 
The central value of the (common) factorization 
and renormalization scale, p,R, for the signal pro- 
cesses has been chosen equal to M H +, whereas 
that of the background processes to 2m t - We 
have checked the uncertainty of the signal due to 
variations of /i/j in the interval M H +/2 < fin < 
2M H +. We find that individual sub-channels 
show a stronger dependence than the total cross- 
section, which varies 5a ~ 28% at the Teva- 
tron II and 5a ~ 18% at the LHC. We have 
also compared our CTEQ4 results with the ones 
obtained with the MRST (central gluon) PDFs [ 
p"H - Again, significant deviations appear for some 
of the individual sub-channels (up to ~ 50%), 
but they are compensated in the sum, leaving a 
5 — 10% thanks to the substraction procedure in 

(§)• 

4. Supersymmetric quantum effects 

In Figs. 0a, b we find the regions of the tan/3 — 
Mh+ plane in which the Tevatron and the LHC 
can find (or exclude) the existence of the charged 
Higgs boson. The curves correspond to the vari- 
ous sets of MSSM parameters indicated in Table 
H For the Tevatron we have included a QCD K- 
factor of 1.5, whereas for the LHC we have just 
set K = 1. From these figures it is obvious that 
the presence of the SUSY corrections alters signif- 
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Figure 1. Discovery regions for the Tevatron and 
the LHC, in the tan/3 — Mjj+ plane, at two levels 
of significance S/ \f~B = 3,5 for the Tevatron, and 
S/\fB = 5 for the LHC. The exclusion limit at 
95% C.L. is also given in each case. We include 
the value of M^o corresponding to each M H + ac- 
cording to the one-loop formulae for the MSSM 
Higgs sector. The tree-level prediction can be 
compared with the SUSY results for the various 
parameter sets of Table @. 



icantly the Higgs boson discovery potential of the 
hadron colliders. It is also patent the Tevatron 
could have a chance to find a charged Higgs bo- 
son in the intermediate range M H ± < 280 GeV. 

Let us now sketch how the various quantum 
corrections have been computed. Among the 
plethora of possible corrections we disregarded 
virtual supersymmetric effects on the gqq and ggg 
vertices and on the gluon propagators. We expect 
those to be of order (a s /47r) • (s/Mg USY ) and 
thus suppressed by a non-enhanced (i.e. tan/3- 
independent) MSSM form factor coming from the 
loop integrals. Therefore, we can neglect these 
contributions as we are only considering effects of 
the form (a s /4ar) n ■ tan™ (3 at large tan/3. The 
cross-section for the signal increases steeply with 
tan/3 and becomes highly significant for tan/3> 
30, while it is much smaller for tan/3 in the low 
interval 2 — 20 where the remaining SUSY correc- 
tions are of the same order or even dominant, so 
our approximation is well justified. Similarly, we 
neglect all those electroweak corrections in ver- 
tices and self-energies which are proportional to 
pure SU(2)l x U(1)y gauge couplings; in partic- 
ular, vertices involving electroweak gauge bosons 
and those involving electroweak gauginos. Fur- 
thermore, we have checked that vertices involv- 
ing Higgs bosons exchange yield a very tiny over- 
all contribution, due to automatic cancellations 
arranged by the underlying supersymmetry. Fi- 
nally, there are the strong gluino-squark diagrams 
and the tan /3-enhanced higgsino-squark vertices 
implicit in chargino-neutralino loops. We have 
extracted the parts of these interactions which are 
(by far) the more relevant ones at high tan (3 and 
confirmed that the remaining contributions are 
negligible. In practice this means that we may 
concentrate our analysis on the interval tan/3 > 
20 where we can be sure that our approximation 
does include the bulk of the MSSM corrections 
while at the same time the cross-section of pro- 
cess (|l]) starts to be sufficiently large to consider it 
as an efficient mechanism for charged Higgs boson 
production. The leading contributions are similar 
to those found in Ref. [ |l2) for the decay process 
t — > H + b, where we also refer for detailed renor- 
malization issues. Recently these dominant ef- 
fects have been conveniently described through an 
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effective Lagrangian approach that contains effec- 
tive couplings absorbing both the leading SUSY 
contributions and the known part of the QCD 
corrections [ [l3| . At high tan j3 the most relevant 
piece is the effective iWJ + -coupling as it carries 
the leading part of the (appropriately resummed) 
quantum effects: 



C 



gVtb rn b (p R ) tan/3 



y/2M } 



w 



1 + Am b 



H+t L b B + h.c. (9) 



The quantity Am b above- see ([12, |13[)- con 



tains the bulk of the supersymmetric contribu- 
tions. Although the leading corrections have been 
identified, we have computed the full set of one- 
loop SUSY diagrams for the relevant tbH + ver- 
tex, which involve similar bunch of diagrams as 
in the on-shell case [ [l^] , but here we have got to 
account for the off-shell external lines, which is a 
non-trivial task. The result is that these off-shell 
corrections amount to a few percent contribution 
that we have included in our numerical analysis. 

Coming back to Figs. |l| and Table [2| let us 
remark that parameter set B constitutes a typi- 
cal case for moderately positive corrections. On 
the other hand set A is over-optimized, and sets 
C and D give more and more negative correc- 
tions. In all cases these corrections decouple very 
slowly with the gluino mass, for fixed values of 
the other masses. Furthermore, they exhibit a 
non-decoupling behavior when all SUSY param- 
eters are scaled up by a common factor [ 12 , |l3) . 
This SUSY non-decoupling behavior [ [l| is associ- 
ated to the properties of the quantity Am b enter- 
ing eq.(||). The violation of the decoupling the- 
orem [ |l4| is caused by the dimensionful Higgs- 
quark-squark coupling, which increases arbitrar- 
ily with the SUSY scale. Through this coupling 
the MSSM Higgs doublet that interacts only with 
the top quark at the tree-level, effectively couples 
to the bottom quark at one-loop. Therefore, that 
interaction cannot be re-absorbed in the param- 
eters of the low-energy Higgs sector after inte- 
grating out the sparticle masses. This welcome 
feature is at the basis of the large radiative cor- 
rections found in our study of the charged Higgs 
production process ([!]), and it could be respon- 
sible for the eventual finding of this scalar bo- 
son. These facts were amply exploited in Ref. [ 



[15| where the suggestion of the present calcula- 
tion was first made and the importance of its cor- 
relation with the neutral Higgs production pro- 
cesses pp(pp) — > hbb + X(h = h°,H°,A°) was 
first emphasized. Such effects and correlations 
could eventually lead to the discovery of SUSY. 
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